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Abstract 

In this paper, a new and efficient algorithm is developed for attitude determination from Global 
Positioning System signals. The new algorithm is derived from a generalized nonlinear predictive filter 
for nonlinear systems. This uses a one time-step ahead approach to propagate a simple kinematics 
model for attitude determination. The advantages of the new algorithm over previously developed 
methods include: it provides optimal attitudes even for coplanar baseline configurations; it guarantees 
convergence even for poor initial conditions; it is a non-iterative algorithm; and it is computationally 
efficient. These advantages clearly make the new algorithm well suited to on-board applications. The 
performance of the new algorithm is tested on a dynamic hardware simulator. Results indicate that the 
new algorithm accurately estimates the attitude of a moving vehicle, and provides robust attitude 
estimates even when other methods, such as a linearized least-squares approach, fail due to poor initial 
starting conditions. 
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Introduction 

Phase difference measurements from Global Positioning System (GPS) signals provides a novel 
approach to determine the attitude of a vehicle. This approach has been successfully applied to air, 1 sea, 2 
and space 3 ' 4 based vehicles. The problem of finding the attitude of a vehicle using GPS signals 
essentially involves a two-step process. First, since phase differences are used, the correct number of 
integer wavelengths between a given pair of antennas must be found. This problem can generally be 
solved using static integer searches or using motion based techniques. Much attention has been placed 
on resolving the integer ambiguity problem over many years (e.g., see Refs. [5-6]). Once the integer 
ambiguities are resolved, then the attitude problem must be solved. The solution to this problem poses a 
difficult task, and has just recently gained attention in the research community. 

The most widely used techniques for attitude determination involve methods that solve Wahba’s 
problem. 7 This problem involves finding a proper orthogonal matrix that minimizes the scalar weighted 
norm-squared residual between sets of 3 x 1 body vector observations and 3 x 1 inertial observations 
mapped into the body frame. Many methods have been developed that solve this problem accurately 
and efficiently (e.g., see Refs. [8-9]). However, the GPS observation is not in the form of a vector 
observation, so finding the attitude using GPS signals is inherently more difficult. 10 

Minimizin g the GPS loss function can be accomplished by using nonlinear least-squares or gradient- 
based search techniques. However, these methods may require a large number of iterations to converge, 
and are not efficient. 11 Cohen’s linearized approach 12 involves finding a small angle rotation which 
maps an initial attitude estimate to the desired attitude matrix. This approach works well for a good 
initial guess, but is not guaranteed to converge to the correct solution for large initial errors. Other 
methods convert the GPS loss function into Wahba’s form. 10,13 The transformation has been shown to be 
exact only when the baselines or sightlines form an orthonormal basis. Significant errors may arise if 
this condition is not true. An extreme example of this scenario is when three baselines are coplanar. 

In this paper, a new and efficient algorithm is derived which determines the attitude using GPS 
observations. The new algorithm is based on a predictive filter scheme for nonlinear systems first 
introduced by Crassidis and Markley. 14 This scheme uses a recursive (one time-step ahead) method to 
“predict” the required model error so that the propagated model produces optimal estimates. The filter 
developed in this paper is essentially reduced to a deterministic approach, since the corrections required 
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to update the model are not weighted in the loss function. The specific name of the new algorithm using 
GPS signals is ALLEGRO (Attitude-Lean-Loping-Estimator, using GPS Recursive Operations). The 
main advantages of the ALLEGRO algorithm over previously developed methods are: 

1) The algorithm is not iterative. 

2) It always converges to the correct solution provided that there is a minimum number of baselines 
and sightlines. 

3) The algorithm is easy to implement. 

An attitude error covariance expression from the general GPS loss function has been developed by 
Crassidis and Markley. 10 It will be shown that the ALLEGRO algorithm produces estimates that have 
exactly the same error covariance provided that the observation sampling is fairly frequent. Therefore, 
the ALLEGRO algorithm minimizes the general GPS loss function. 

The organization of this paper proceeds as follows. First, the concept of the GPS phase difference 
observation is introduced. Then, the general loss function used for GPS-attitude determination is 
reviewed. Next, for completeness the optimal attitude error covariance derivation is shown. Then, the 
generalized predictive filter for nonlinear systems is reviewed, followed by an application of this scheme 
to the GPS loss function. Also, an attitude error covariance expression is derived for the predictive 
attitude determination algorithm. Finally, the algorithm is tested using a GPS hardware simulator. 

Background 

In this section, a brief background of the GPS phase difference measurement is shown. The main 
measurement used for attitude determination is the phase difference of the GPS signal received from two 
antennas separated by a baseline. The wavefront angle and wavelength are used to develop a phase 
difference, as shown in Figure 1 . The phase difference measurement is obtained by 

b i cos 0 = X(Atf>-n) (1) 

where b t is the baseline length (in cm), 0 is the angle between the baseline and the line of sight to the 
GPS spacecraft, n is the number of integer wavelengths between two receivers, A<p is the phase 
difference (in cycles), and X is the wavelength (in cm) of the GPS signal. The two GPS frequency 
carriers are LI at 1575.42 MHz and L2 at 1227.6 MHz. As of this writing, non-military applications 
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generally use the LI frequency. Then, assuming no integer offset, the phase difference b.<f> can be 
expressed by 

&<t> = b T As (2) 

where s elR 3 is the normalized line of sight vector to the GPS spacecraft in an inertial frame, b eiR 3 is 
the baseline vector in wavelengths, which is the relative position vector from one antenna to another, and 

A is the attitude matrix, which is an orthogonal matrices with determinant 1 (i.e., A T A = I). 

Attitude determination using GPS signals involves finding the proper orthogonal matrix A that 
minimizes the following generalized loss function 


. m n 

( 3 ) 

»•=! J = 1 

where m represents the number of baselines, n now represents the number of observed GPS spacecraft, 
A <f> denotes the phase difference measurement, and ay denotes the standard deviation of the y ,th 

measurement error, which is assumed to be a zero-mean stationary Gaussian process. The standard 
deviation is 0.5cm/U = 0.026 wavelengths for typical phase noise. 12 

A convenient parameterization of the attitude matrix is the quaternion representation, defined as 15 

<7 = 1 ^ 13 ] (4) 

“ L*4j 

with 


?13 " 


?3 


= |sin(0/2) 


(5a) 


? 4 =cos(0/2) (5b) 

where e is a unit vector corresponding to the axis of rotation and is the angle of rotation. The 
quaternion satisfies a single constraint, given by 
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The attitude matrix is related to the quaternion by 



(7) 


with 



<7473x3+ [^x] 




+[2,3 *] 



(8a) 


(8b) 


where 73*3 is a 3 x 3 identity matrix. The 3x3 matrix X J is referred to as cross product matrix 
since axft = [a x]6, with 


[gx]- 


0 

a 3 


~ a 3 a 2 
0 -a. 


-a 2 a, 0 


(9) 


From Equation (3) it is clear that the quaternion representation leads to a loss function that is quartic in 
the quaternions. This is not equivalent to the familiar attitude determination loss function posed by 
Wahba. 7 In fact, a conversion to Wahba’s problem has been shown be optimal only when the baselines 
or sightlines form an orthonormal basis. 10 Therefore, in general, the GPS loss function poses a more 
difficult problem to solve than the standard vector-observation loss function in Wahba’s problem. 

An attitude error covariance can also derived from the GPS loss function in Equation (3). This is 
accomplished by using results from maximum likelihood estimation. 10,16 The Fisher information matrix 
for a parameter vector x is given by 


F XX =E\ 


dxdx k 


^(*) 


( 10 ) 


where E{ } denotes expectation, and J(x) is the negative log likelihood function, which is the loss 
function in this case. If the measurements are Gaussian and linear in the parameter vector, then the error 
covariance is given by 
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( 11 ) 


p - F~ 

1 XX 1 XX 

Now, the attitude matrix is approximated by 

A = e~t— X ^ true *(/ 3x3 -[^ x ] + |[^a x f]^true ( 12 ) 

where 5a represents a small angle error (for the quaternion 25^^^ 5a). Equation (12) is next 

substituted into Equation (3) to determine the Fisher information matrix. First-order terms vanish in the 
partials, and third-order terms are small because we assume the probability distribution to be 
approximately symmetric about the mean. Also, asstiming that the quartic terms are negligible (see [17] 
for a Gaussian approximation to fourth-order terms) leads to the following form for the optimal 
covariance 

m n 

f o P t * 

,WJ=I 

Note that the optimal covariance requires knowledge of the attitude matrix. However, if the baselines 
are non-coplanar then the optimal covariance can be determined without the attitude knowledge. 10 

There are a number of methods available to minimize the GPS loss function shown in Equation (3), 
including the standard parameter optimization techniques, such as the gradient method. 18 However, 
these methods are usually computationally inefficient. A more practical approach uses a linearized least- 
squares method. 12 This begins by performing a first-order linearization about a nominal attitude, so that 

A = ^o(/3 X 3 +[£#*]) (14) 

where Aq represents some nominal attitude, and 59 represents a small angle correction. Then, defining 
a perturbation equation for the phase difference measurement leads to 

5<p- A^q = s T Aq[b x]59 + v = h T 59 + v (15) 

where v represents the Gaussian measurement noise. Equation (15) represents a linearized sensitivity 
equation between the measured differential carrier phase and the perturbation to the initial attitude guess. 
All available differential phase measurements can be stacked into a single linearized vector equation, 
given by 
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(16a) 


'hi i 

H = : 

—mn 

50 = H56+v (16b) 

Therefore, Equation (16b) can be used to find a least-squares estimate of the attitude from the nominal 
attitude. In practice, the solution is usually obtained iteratively by using the previous epoch as an initial 

guess. Also, it is easy to see that ( H T R~ l H ) is equivalent to the attitude error covariance expression 

in Equation (13), where R is the diagonal covariance matrix of the measurement error process y. The 
linearized approach provides an efficient method for attitude determination; however, it is sensitive to 
the initial attitude guess, which may cause divergence problems (as will be shown). 

Predictive Attitude Determination 

In this section, the ALLEGRO algorithm is derived using a nonlinear predictive approach. First, a 
brief review of the nonlinear predictive filter is shown (see Ref. [14] for more details). Then, the filter 
algorithm is reduced to a deterministic-type approach for attitude determination. Finally, a covariance 
expression for the attitude errors using the ALLEGRO algorithm is derived. 

Predictive Filtering 

In the nonlinear predictive filter it is assumed that the state and output estimates are given by a 
preliminary model and a to-be-determined model error vector, given by 

x{t) = f{x{t), t) + G(t)d(t) (17a) 

y(t) = £(i(0>') (17b) 

where / is the model vector, x(t) e 91 p is the state estimate vector, d(t) €*R ; is the model error 

vector, G(t) e yi px! is the model-error distribution matrix, c e is the measurement vector, and 
y(t) e is the estimated output vector. State-observable discrete measurements are assumed for 
Equation (17b) in the following form 

K^) = eU(^).^)+i:(r yfe ) (18) 
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where y(/*) e'-R 7 ” is the measurement vector at time t k , x(t k )e9i p is the true state vector, and 

v(t k ) is the measurement noise vector which is assumed to be a zero-mean, stationary, Gaussian 
white-noise distributed process with 

£{v('*)} = 0 (19a) 

E {l{tk h T (^ )} = R s kk ( 1 9b) 

where R e y\ mxm is a positive-definite covariance matrix. 

A loss functional consisting of the weighted sum square of the measurement-minus-estimate 
residuals plus the weighted sum square of the model correction term is minimized, given by 

• / = |{i('*+l)-^'*+l)} r *~ 1 + d{l t ) (20) 

where W e is weighting matrix. The necessary conditions for the minimization of Equation (20) 
lead to the following model error solution 

*t) = -{[A(A»)s(i i )] 7 'jr 1 A(A»)s(i i ) + ir}“ , [A(A/)s(i t )] r «- 1 [ 2 (£ t ,A/)-y(, t+l ) + ^(, i )] (2 i) 

where At is the measurement sampling interval, S(x) el H mxl is a generalized sensitivity 

matrix, and A(A/) € <R mxm is diagonal matrix with elements given by 


A it - - 


A t Pl 


Pi 


.1 


i = 1,2 


( 22 ) 


where p h i = l,2,...,m, is the lowest order of the derivative of c, (x(t)) in which any component of d(t) 

first appears due to successive differentiation and substitution for £,•(/) on the right side. The /' th 
component of z(x, At) is given by 


z,(x,At) = ^^-L k / (c i ) 


*=1 


where £y-(c,) is the £ th Lie derivative, defined by 


(23) 
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L g,[ L f V/)]}’ i = h2,...,m 


(25) 


The I th 


row of 5(i) is given by 


s, = 



where is the y ,th column of G(t), and the Lie derivative is defined by 

r . t 3L p r\ct) 

L zA L f ^r df g J' j = ( 26 ) 

Equation (26) is in essence a generalized sensitivity matrix for nonlinear systems. Therefore, given a 
state estimate at time t k , then Equation (21) is used to process the measurement at time to find the 
d(tfr) to be used in to propagate the state estimate to time t k+i . The weighting matrix W serves 

to weight the relative importance between the propagated model and measured quantities. If this matrix 
is set to zero, then no weight is placed on minimizing the model corrections so that a memoryless 
estimator is given. 

ALLEGRO Algorithm 

In the ALLEGRO algorithm it is assumed that the model is given by the quaternion kinematics 
model. This algorithm requires no dynamics model; it assumes that the attitude rate is adequately 
modeled by a constant model error d between measurements, so that 

l = (27) 

where q denotes the determined quaternion. Since the phase difference measurements are used as the 
required tracking trajectories, the output vector in Equation (18) is given by (dropping the subscript ij 
for the moment) 

£(i) = b T A[q^s (28) 

The lowest order time derivative of q in Equation (28) in which any component of d first appears in 
Equation (27) is one, so that Pj = \. For a deterministic attitude solution (i.e., a memoryless approach) 



the weighting matrix \ V is set to zero in Equation (21). The remaining quantities in Equation (21) can 
be shown to be given by 


A = At /3 X 3 


kj[A{q)s_f x] 



X 

1 

bl[A(q)s% x] 



U*ll 

A 

y = 

hi 


- 

h T m A[q)sJi 

_T 


T 

Z('* + i) = [A^ A 1, •••, M>mn] 


R = diagftr i ^ .... <r mn ] 
A/) = 0 


(29a) 

(29b) 

(29c) 

(29d) 

(29e) 

(29f) 


where the superscript A denotes that the quantity is measured at time /* +1 (all other quantities in these 
equations are at time fy). Therefore, the following model error equation is developed 


i-l 


d{t k >dk = - 


A/l 


m n 


*\ b -i£\ As -j *] 

;=1 j-l 


m n 


x M a $ < 30 > 

i=l 7=1 


It should be noted that Equation (30) represents an exact linearization for an interval A/. 19 However, for 
practical applications the sampling interval should be well below Nyquist’s limit. 20 The determined 
quaternion can be found by integrating Equation (27) from time to /^ + j. Since d is assumed constant 
over this interval, a discrete propagation for Equation (27) can be used, given by 




where 



(31) 


(32a) 


to 



(32b) 


fik s cos^|^||A/J 
Yk 

P k =d k l\\d k \\ 


(32c) 

(32d) 


In order to derive an attitude error covariance from Equation (27), a propagated expression can be 
derived using a similar approach found in Ref. [21]. The attitude error equation is given by 

5a = -[d x\Sa + Sd (33) 

where 5d is a model error perturbation. The discrete propagation is given by 

= ®k dpLk + —k (34) 

where 


(35a) 

r i= f 4 V [ - X 1'<* (35b) 

JO 


Next, the true output is given by using a first-order expansion in the predictive filter output, 14 so that 


yn+i-yjc + AlSkdjt + Vk+l ( 36 > 

where S and d_ correspond to true quantities of S from Equation (29b) and d from Equation (21), 
respectively. Therefore, the model error is given by 



(37) 


where 


K k =(slR-'s k )~' Sir' 

Next, using a small angle perturbation in the attitude matrix, similar to Equation (12), leads to 

y_ k -y k ~s k 5a k 


(38) 


(39) 


it 



Now if Sa k is small, using the right-hand side of Equation (29b), the following approximation can be 
used 

S k » S*(/ 3x3 + [§a k x]) (40) 

Therefore, since K k S k = / 3 X 3 ,the model error equation is now given by 

dk ~^ik/ At + K ky.k+ l/ A/ + ( / 3x3 + (41) 

Using the fact that Sd k -d_ k —d_ k leads to the following error angle equation 

te k+l =<t>kfc_k- T k ScLk!^ - r * K k Yk+i/to + r* [Ik x]5a k (42) 

If At is small, as assumed in this approach (i.e., the sampling interval is well within Nyquist’s limit, 
Ifc I* < nf 10), 20 then the quantities in Equation (35) can be approximated adequately by 

0^(7 3x3 -A/[^ x]) (43a) 

tl 3x3 (43b) 

Substituting these quantities into Equation (42) leads to 

dSLk + i"-K k Yk + i (44) 

The cancellation of the terms in Sa k reflects the fact that setting W =0 in Equation (21) gives a 
memoryless estimator. Now the attitude error covariance is given by 

p k + 1 s E {&Lk + 1 &!+ 1} = K k RK I (45) 

Therefore, from the definitions of S k , K k , and R, the attitude error covariance expression for the 
ALLEGRO algorithm becomes 

m n 

= YL a t 

./= 1 7=1 

Note that the attitude matrix in Equation (46) is evaluated at time t k and that the sightlines are given at 
time t k+ This may be simplified by using the following attitude propagation which is valid for small 
At 


[ A (ii) s -j x ]htf[ A {ik) s -j X ] 1 
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(47) 


A {ik + |) = ( 7 3x3-Afe 

The inverse recursion for /f^ j can be adequately approximated by 

A {ik) * ( 7 3 x 3 + x ] ) A {ik + \) < 48 > 

Substituting Equation (48) into Equation (46) leads to 

r i-l 

m n j- 

YL a t{tj\Wi 

.w>i 

where 

4 s [& + to[dt x])x (* +1 4 (50) 

The term in Equation (49) that involves [d k x] is typically three orders of magnitude less than the term 
that doesn’t involve \d k x], and the term that is quadratic in [d k x] is typically six orders of magnitude 
less than the term that doesn’t involve \d k x] . Thus, Equation (49) reduces down to 

m n 

^+i* 

./= 1 7=1 

Therefore, the attitude error covariance at time t k is given by 

r t-1 

m n T 

p~ YL°f\ A ij*]^ As -A (52) 

./=! 7=1 

This expression is equivalent to the optimal covariance shown by Equation (13). Therefore, the 
ALLEGRO algorithm is in essence equivalent to solving the generalized loss function in Equation (3). 
Although the approximation in Equation (52) is valid only for small A/, this poses no problem for 
typical on-board applications (e.g., for a typical vehicle in low-Earth orbit undergoing motion of one 
revolution-per-orbit, a sampling interval of 100 seconds is more than sufficient for Equation (52) to be a 
valid approximation). Also, the inverse in Equation (30) is sufficient to determine P k +\, as shown by 
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Equation (46)-(51). Therefore, the ALLEGRO algorithm inherently computes the attitude error 
covariance as part of its solution. Finally, Ref. [22] shows an analysis of robustness with respect to 
initial condition errors. It is shown that the estimated error in predictive filter is always bounded for any 
initial condition, which makes the ALLEGRO more robust than a linearized least-squares algorithm. 

There are many advantages of the ALLEGRO algorithm over previous methods. These advantages 
include: 

1) The ALLEGRO algorithm can provide estimates even when the baselines are coplanar, which is 
an advantage over the methods shown in Refs. [10] and [13] that convert the GPS problem into a 
form equivalent to Wahba’s problem. Also, it has been shown in Ref. [10] that the attitude of a 
vehicle can be determined with a minimum of two baselines and two sightlines (to within a sign 
change). This is also true for the ALLEGRO algorithm, for which the solution will converge to 
the true attitude as long as the initial condition is in the correct hemisphere. 

2) Unlike gradient based-methods the ALLEGRO algorithm is non-iterative, which provides a more 
numerically stable algorithm. 

3) The ALLEGRO algorithm is robust with respect to initial condition errors, which is an advantage 
over the linearized least-squares algorithm. 

4) The computational burden of the ALLEGRO algorithm is low, since the algorithm is easily 
programmable using Equations (30) and (31). 

Hardware Simulation 

A hardware simulation of a typical spacecraft attitude determination application was undertaken to 
demonstrate the performance of the ALLEGRO algorithm. For this simulation, a Northern Telecom 40 
channel, 4 RF output STR 2760 unit was used to generate the GPS signals that would be received at a 
user specified location and velocity. The signals are then provided directly (i.e., they are not actually 
radiated) to a GPS receiver that has been equipped with software tracking algorithms that allow it 
operate in space (see Figure 2). 

The receiver that was used was a Trimble TANS Vector; which is a 6 channel, 4 RF input 
multiplexing receiver that performs 3-axis attitude determination using GPS carrier phase and line of 
sight measurements. This receiver was modified in software at Stanford University and NASA’s GSFC 
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to allow it to operate in space. This receiver model has flown and operated successfully on several 
spacecraft, including: REX-II, OAST-Flyer, GANE, Orbcomm, Microlab, and others. 

The simulated motion profile was for an actual spacecraft, the Small Satellite Technology Initiative 
(SSTI) Lewis satellite, which was launched on August 22, 1997 (see Figure 3). The orbit parameters 
used for the simulation are given in Table 1. This mission actually carried a GPS attitude determination 
experiment to assess the performance of the GPS attitude measurements on-orbit. Although the 
spacecraft was lost due to a malfunction not related to the GPS experiment shortly after launch, this 
motion profile is nonetheless very representative of the types of attitude determination applications that 
are found on satellites. 


Table 1. SSTI Lewis Orbit parameters 


Semimajor axis (a) 

6901.137 km 

Inclination (i) 

97.45 deg 

Right Ascension of Ascending Node (RAAN) 

-157.1 deg 

Eccentricity (e) 

0.0001 

Pointing profile 

Earth pointed 

Launch date 

August 22, 1 997 


The antenna separation distances are 0.61 m, 1.12 m, and 1.07 m, respectively. One antenna (in 
baseline 3) is located 0.23 m out of plane (below) the other three antennas. On the spacecraft, the 
antennas are mounted on pedestals with ground planes to minimize signal reflections and multipath. For 
the simulation, the signal was provided to the GPS receiver without multipath noise. The simulated 
SSTI Lewis spacecraft has four GPS antennas that form three baselines. The baseline vector 
components in wavelengths are given by 



‘ 2.75 ‘ 


" 0.00 ■ 


'-3.93' 

*1 = 

1.64 

-0.12 

> —2 “ 

6.28 

-0.17 

> —3 = 

3.93 

-1.23 


(53) 


Quantities such as line biases and integer ambiguities are first determined before the attitude 
determination algorithms are tested. The GPS raw measurements are processed at 1 Hz over a 40 minute 
simulation. A plot of the number of available GPS spacecraft for the simulated run is shown in Figure 4. 
During the beginning of the run there are 5 to 6 available spacecraft. At the end of the simulation this 
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drops down to about 4, which means that a degraded performance is possible (this also depends on the 
geometry of the spacecraft, see Ref. [12] for Geometric Dilution of Precision). 

For the first simulation the sightlines and baselines are used to form simulated phase difference 
measurements with Gaussian measurement errors. This is not a totally realistic simulation; however, it 
is useful to quantify the effectiveness of the ALLEGRO algorithm. A plot of the (roll, pitch, yaw) 
attitude errors with 3-sigma outliners using Equation (46) is shown in Figure 5. Clearly, the ALLEGRO 
algorithm provides estimates that agree with the optimal standard deviation predictions. 

The remaining runs use the actual phase measurements from the receiver. This provides a more 
realistic scenario. The linearized least-squares approach using Equations (14)-(16) is also used to 
determine the attitude. A plot of the determined attitude using the ALLEGRO and least-squares 
algorithms is shown in Figure 6 (the glitch between 10 and 15 minutes is due to receiver outages). The 
ALLEGRO attitudes exactly match the least-squares determined attitudes. In order to test the robustness 
of the both algorithms, each is started with a poor initial attitude guess. A plot of the attitude errors 
during the iteration stage of the least-squares algorithm is shown in Figure 7. Clearly, the least-squares 
algorithm does not converge to the correct solution. This is due to the small angle approximation in 
Equation (14). The same initial condition is applied to the ALLEGRO algorithm. Since the ALLEGRO 
is sequential and non-iterative, convergence is given over sampled intervals. A plot of the attitude errors 
is shown in Figure 8. Clearly, the ALLEGRO algorithm converges to the correct solution (after 3 
sampling intervals for this initial condition). To further test the robustness of the ALLEGRO algorithm, 
a Monte Carlo analysis has been performed using 1000 normalized random initial conditions. A plot of 
the convergence rates is shown in Figure 9. In all cases, convergence is achieved within 19 sampling 
intervals (most converge within about 10 sampling intervals). 

Finally, a test has been performed on the computational efficiency of the ALLEGRO algorithm. The 
number of floating point operations (FLOPS) has been evaluated using MATLAB. Both methods 
calculate the attitude error covariance as part of their solutions. A comparison with the least-squares 
algorithm is slightly misleading, since ALLEGRO is non-iterative. It has been determined that the only 
major difference between them is the ALLEGRO algorithm propagates a quaternion model. However, 
the computational expense of this propagation is smaller than 75 FLOPS, which is almost an order of 
magnitude less than doing a second iteration in the least-squares algorithm (even for only the two 
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baseline and two sightline case). Therefore, the ALLEGRO algorithm is computationally comparable or 
better than the least-squares algorithm. 

Conclusions 

In this paper, a new optimal and efficient algorithm has been developed for attitude determination 
using Global Positioning System signals. It has been shown that the standard GPS loss function is 
inherently difficult to solve. The new non-iterative algorithm provides sequential estimates using a 
recursive one-time step ahead approach. Attitude determination is accomplished by determining the 
angular velocity components used to propagate a simple quaternion kinematics model. An attitude error 
covariance expression has been derived for the new algorithm. This covariance has been shown to be 
equivalent to the optimal covariance, derived from maximum likelihood, if the sample interval is small 
enough (which poses no problem for most applications). The algorithm was tested on a hardware 
simulator using an actual receiver. Results indicated that the new algorithm is computationally 
comparable to a linearized least-squares approach, while providing robustness with respect to initial 
conditions error. Therefore, the algorithm is exceptionally suitable for on-board applications. 


Appendix: Alternative Covariance Derivation 

In this section another approach for the attitude error covariance in the ALLEGRO algorithm is 
derived. Linearizing Equation (31) in d k gives 



(Al) 


where the identity Cl(d k )q k =^q^d k was used. Next, use A = (/ 3 X 3 - [ 5a x]) , so that 

bj Asf = bj. .<4 true sf - bJ[Sa xj^rue sf 

= &0 i j-Sa T ri i j 


(A2) 


where = j^Q-ue sf and A A/ denotes the true phase difference. Using A Ay - A Ay =v < y, and 
substituting Equations (A2) and (30) into Equation (Al) yields 


1**1 =[/4*4 - jKs*) 


1-1 


m 


TL 

m y=i 


-2 T 
cr. : r] .7 7 '. 

U ±ij±ij 


m n 

<7 


2.2.' 

/=i y=i 


V U-H V ij 


(A3) 
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Now, using the fact that the true quaternion ((7) can be represented by 

£ = /4 x 4~ q (^)]? (A4) 

yields the following approximation to within first-order 

r n-l 

m n m n 

=£*-2^*) ZZ^Ay ZZ^AA (A5> 

_ 1=1 7= 1 _ 1=1 y=l 

Defining Sq = q -q gives the following quaternion error covariance 

r l-l 

m n 

ZZ 'tf’hjll zT (h ) 

/=1 y=l 

Therefore, using the same principles for the attitude error covariance derivation in Ref. [8] and from the 
analogy in Equations (46) through (51) gives 

m n 

YL°t 

/=1 7=1 

which is again the same expression as in Equations (13) and (52). 
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